Single-particle Excitation Spectra of Ceo Molecules and Monolayers 
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In this paper we present calculations of single-particle excitation spectra of neutral and three- 
electron-doped Hubbard Ceo molecules and monolayers from large-scale quantum Monte Carlo sim- 
ulations and cluster perturbation theory. By a comparison to experimental photoemission, inverse 
photoemission, and angle-resolved photoemission data, we estimate the intermolecular hopping in- 
tegrals and the Ceo molecular orientation angle, finding agreement with recent X-ray photoelectron 
diffraction (XPD) experiments. Our results demonstrate that a simple effective Hubbard model, 
with intermediate coupling, U = 4t, provides a reasonable basis for modeling the properties of Ceo 
compounds. 
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I. INTRODUCTION 

Cgo compounds occupy a unique place within the pan- 
theon of high-temperature superconductorsi 2 -^^ One 
would not normally expect such high transition temper- 
atures (30-40K) to arise from the conventional phonon 
mechanism of superconductivity. However, the unusual 
molecular structure of these compounds results in higher 
frequency phonons and stronger couplings than are pos- 
sible in typical elemental compounds and alloys.— In a 
similar way, the effects of electron-electron interactions 
on these molecules are more complex and, it has been 
argued^ could also provide the mechanism for attrac- 
tive pairing. There is general agreement that the attrac- 
tive mechanism is intramolecular, i.e., it arises from in- 
tramolecular vibrations or from Coulomb correlations on 
a molecule or, perhaps most likely, from a combination 
of the two effects.— This intramolecular attraction, to- 
gether with the narrow bands and correspondingly large 
density of states (DOS) for intermolecular hopping, is 
believed to account for the observed high transition tem- 
peratures. The two contrasting theoretical approaches, 
which emphasize either phonons or Coulomb interac- 
tions, can be roughly characterized as an LDA electronic 
structure coupled to Jahn- Teller phonons, compared to 
a Hubbard model of the Ceo molecule with tight-binding 
intermolecular hopping. 

Recently Yang et alM- reported angle-resolved pho- 
toemission (ARPES) and LDA calculations of the band 
structure of two-dimensional (2D) hexagonal monolayers 
of Ceo and K3C60 on the surface of silver, and Wehrli 
et al*£ have proposed an electron-phonon mechanism to 
explain the DOS observed in Yang et al.'s study as well as 
spectra from earlier experiments . 13 i 14 i 15 Since their cal- 
culations do not capture the distinctive correlation effects 
of the Hubbard model within a molecule, it is of inter- 
est to examine how the DOS and other properties are 
affected by such correlations. 

In this paper, we calculate the photoemission spec- 
tra (PES) and inverse photoemission spectra (IPES) 
of 2D hexagonal Ceo monolayers including the strong, 




FIG. 1: (Color online) 2D hexagonal monolayer of identically 
oriented Ceo molecules (lattice constant a = 10.02A for the 
superlattice 1 !) and its corresponding first Brillouin zone (BZ). 
The definition of the molecular rotation angle 9 is the same 
as in the caption of Fig. 4 of Ref. [Til , which is defined by two 
lines, both projected onto the plane of the layer, one from a 
molecular center to a NN molecular center, as shown, and the 
other from a molecular center through a pentagon face cen- 
ter. 8 — 0° when these two lines coincide. Counterclockwise 
rotation of the molecule corresponds to a positive rotation an- 
gle. For K3C60, the black squares represent two K + ions, one 
above the other, while gray ones denote a single K + ion.— 



atomic, on-site electron-electron Coulomb interaction U 
within a single molecule. The calculation is performed 
by cluster perturbation theory (CPT ) 6 i 17 using quan- 
tum Monte Carlo (QMC) data3 This method, which 
we call QMCPT, is ideally suited to calculating the elec- 
tronic properties of solids composed of complex units 



such as fullerene molecules. To distinguish, we shall refer 
to the standard CPT approacbJ&il as EDCPT since it 
uses exact diagonalization (ED) results. By comparing 
the calculated single-particle excitation spectra, for inter- 
mediate interaction strength, U = At, with experimental 
data t 11 ^ 19 we can estimate the intermolecular hopping in- 
tegrals t and the molecular orientation angle, 9, which 
are useful for an effective Hubbard model Hamiltonian 
description of Cgo compounds. From the best fit param- 
eters, we find a rotation angle of 9 = 64° (as defined in 
Fig. [TJ and assuming that the C 60 molecules are all ori- 
ented identically) which is consistent with experimental 
measurements j 11 ' 20 



II. MODEL 

The Hamiltonian for the 2D hexagonal Cgo superlattice 
is given by 

H = H a + V, (1) 

#o = £*To. (2) 
I 

where 

Hq=~ M c L c j<x + h.c.) + ™;T n a ( 3 ) 

(Iijj)a iel 

is the Hubbard model on the Ith Cqo molecule, and 

v = - E ;/ '"<••'; • • (4) 

(Ii,Jj)a 

is the hopping Hamiltonian between a pair of carbon 
(C) atoms i and j on two nearest neighbor (NN) Ceo 
molecules I and J, respectively. The operator notations 
are standard for the usual Hubbard model. Inside the 
molecule U is the on-site Coulomb interaction energy for 
two electrons on the same C atom, and Uj is the hopping 
integral between two NN carbon atoms. A molecular 
U, parameterizing the interaction energy for two elec- 
trons on the same Cgo molecule, and hence distinct from 
our atomic U, is also discussed in the literature. See for 
instance Ref. [l^. For the Ceo molecule, there are two 
kinds of bonds between NN sites: — t for single bonds 
(1.461), and U 3 = 1.2* for double bonds (1.40l). We set 
t = 2.72eV according to Ref. l2ll. The value of the on-site 
Coulomb energy, U, is taken to be U = 4i, an interme- 
diate coupling value for which the sign problem for the 
single-molecule QMC is manageable* 22 

The shortest distance between two Ceo molecules is 
about 3.0 A. The relative sizes of the direct hopping inte- 
grals Oin y jj between two C atoms on two NN molecules 
are given by 2 ^ 

an.j 3 = [y CT (d)-^(d)](R /r d)(R J ,-d)+V;(d)(R /2 -R Jj ), 

(5) 
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FIG. 2: (Color online) DOS of (a) Ceo and (b) single 
molecules from QMC and MEM for U = 4t and fit = fO. A 
total of 465 and 620 bins of space-time Green's functions are 
collected for Ceo and C|o~, respectively, for MEM analysis. 
Each bin is an average over f 00 space-time Green's functions, 
which are collected after each QMC sweep over the whole 
space-time lattice. 

where R is a unit vector along radial direction, d = d/d 
is a unit vector pointing from atom Ii to Jj, and 

V a {d) = -W n (d) = 4- exp[-(d - da) I L\, (6) 

with L = 0.5051, and d = 3.0l.— We leave the overall 
prefactor t as a parameter to be determined later. 

We see from Fig. [T]that for K3C60 the Cgo molecules 
are separated by K + ions, either two ions, one above the 
other or single ions. The indirect hopping between Cgo 
molecules via K + ions is thus of comparable importance 
to that of the direct C-C hopping integrals. Expressions 
for the indirect hopping between C atoms via the K + 
ions are given by 2 ^ 

.ind _ ^Ji,7^Jj,7 f7 \ 

^ ec — £K 

7 

where ec — &k = — 4eV is the energy difference between 
a C atom and K + ion, 7 = 1, • • • , 9 is K + index, and the 
C-K hopping integral is given by 

t^^lMD^e-^-*™^, (8) 
mcr 

where d is C-K separation, m is electron mass, i? mm 
is the shortest C-K separation for a given K + and Cgo 
molecule, and D — 0.473^ We take the indirect hopping 
matrix elements to be fixed while the direct hopping ma- 
trix elements scale with the fitting parameter t . For the 
value of t used to fit the ARPES data, we find that the 
largest direct hopping matrix elements are substantially 
larger than the indirect hopping matrix elements. 
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We note that our model Hamiltonian, given by Eqs. 
©-(HI), neglects both the long-range Coulomb interac- 
tions and the effect of the silver substrate, both of which 
would be difficult to include in QMC calculations. One 
would expect the long-range Coloumb interaction to in- 
crease the effective molecular U. On the other hand, the 
silver substrate, as well as the other molecules, will screen 
this interaction and, consequently, these two effects may 
partially cancel. In any case, we will see below that the 
simple Hubbard model and parameters employed here 
yield results in generally good agreement with experi- 
ment. 



QMCPT, 2D, OBC 
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III. RESULTS 

A. Single Molecule DOS 

We begin by considering results for a single Cgo 
molecule. QMC simulations were used to calculate 
the imaginary-time Green's functions for Ceo and C|^~ 
molecules, based on the single-molecule Hamiltonian, Hq 
and standard QM C 25 ' 26 ' 27 for given values of t and U. 
The QMC step for the single molecule is by far the most 
time-consuming part of the calculation. The maximum 
entropy method (MEM) was then used to analytically 
continue the imaginary-time Green's functions to real 
frequency^ which gives the single molecule DOS shown 
in Fig. [2] for T = O.lt. From the figure, we see that neu- 
tral Ceo is an insulator with an energy gap about 1.05i or 
2.86eV. The presence of this gap is not surprising since 
neutral Cqq has a filled shell with a sizable excitation gap 
even for U = 0. In fact the gap obtained by extended 
Huckel calculations for the neutral model is also about 
It, although the correlations for that case are very dif- 
ferent than for U — At. For the Cqq molecule, the Fermi 
energy lies in a region of non-zero density of states. 
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FIG. 3: Single particle spectral functions of a 2D Hubbard 
model from (a) QMCPT and (b) EDCPT with open bound- 
ary conditions. The cluster is of dimension 3x4. The two 
methods predict very similar single-particle excitation ener- 
gies and energy gaps. Note that special k points F, M, X are 
for a 2D square lattice, different from those in Fig. [1] 



C. Monolayer DOS 



B. QMCPT for 2D Square Lattice 

The remaining perturbative part of the QMCPT calcu- 
lation is, by comparison, very economical. It uses the sin- 
gle molecule Green's functions calculated above to con- 
struct any superlattice Green's functions, incorporating 
the effect of the intermolecular hopping Hamiltonian, V , 
and different molecular orientations with negligible addi- 
tional computational effort. This method has been ap- 
plied to the simpler case of the 2D Hubbard model on 
a square lattice, which is metallic for U = but which 
develops an insulating Mott-Hubbard gap for sufficiently 
large U. The "cluster" in this case is a small unit of the 
lattice which can be treated either by QMC or by exact 
diagonalization. Results for U — 8t, obtained by QM- 
CPT and by EDCPT, based on a 3 x 4 cluster of sites^ 
are shown in Fig. [3J 



Next we consider the case of Cgo monolayers. For these 
calculations we use values of t and the orientation an- 
gle 9 which are obtained by fitting photoemission and 
diffraction data as will be discussed below. We show the 
DOS for monolayers of both Ceo and K3C60, in panels (b) 
and (c) of Fig. 3J Also shown for comparison, in panel 
(a) of Fig. 21 is the DOS for the same values of t and 9, 
but for U = 0. The lattice of neutral Ceo molecules is 
very different from the 2D square lattice of atomic sites, 
for which spectra are shown in Fig. [31 because the Ceo 
molecule has an excitation gap as was shown in the top 
panel of Fig. H For U = [(a) panel of Fig. g], the C 60 
monolayer is a band insulator with a gap of about O.bt 
which is about half the value for the single-molecule case 
because the bands above and below the gap are brodened 
by temperature effects and the intermolecular Hamilto- 
nian V, thus reducing the gap. We see, from the (b) 
panel where U = At, that the effect of U is to enlarge the 
gap back to a value (~ 2.8eV) slightly greater than It 
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FIG. 4: (Color online) DOS of hexagonal (a) Ceo monolayer 
with U — 0, (b) Ceo monolayer with U = At, and (c) K3C60 
monolayer with U — At from QMCPT calculations for t — 
— 0.3i and T — O.lt. Shaded areas are occupied by electrons. 



due to short-range, Mott-Hubbard corelations. This gap 
value is comparable to the value of 2.3eV measured by 
Lof et al. in Ref. [IjJ. Our calculated gap value may be 
somewhat larger than the experimental gap due to our 
neglect of the silver substrate, which may reduce the en- 
ergy gapi22 By contrast, for the case of K3C60 [(c) panel] 
the monolayer remains metallic with its Fermi energy ly- 
ing within a narrow conduction band (Fig.f?]). 

To determine the parameter t and orientation an- 
gle 9, we vary t and 9 values in the QMCPT calcula- 
tions. The resulting energy-momentum dispersion curves 
(along the YM direction) for the hexagonal monolayer of 
K 3 Cgo are then compared with the experimental one from 
ARPESJi We find that the best fit is obtained when 
t = -0.3* with a rotation angle 9 = 64° ± 4°. From 
Fig. [5J we see that both the experimental and theoretical 
band widths are about lOOmeV. However the theoretical 
band dispersion in the k range [—0.8, — 0.6]A _1 is rather 
weak around the Fermi energy, while in the experimental 
dispersion figure, the band dispersion in that momentum 
range is much more obvious. We think that this results 
from two causes: First MEM has difficulty reconstruct- 
ing the weak features of the band dispersion curve, and, 
second, the experimental scanning direction is not ex- 
actly along TMiii The above determined rotation angle 
9 = 64° ± 4° is in good agreement with XPD structural 
data3 

The densities of states in the lower two panels of Fig. 
4 arc directly comparable to the PES data (Fig. IE in 
Ref. [Tlh of Yang et ai, we draw only the photoemission 
(PE) part of the calculated DOS together with the ex- 
perimental PES data in Fig. [6] We see that the QMCPT 
results agree reasonably well with the PES data close 
to the Fermi energy, with approximately the same en- 
ergy gaps for the Cgo monolayer and similar peaks at the 
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FIG. 5: Determination of NN molecular hopping integral t 
and molecular orientation angle 9 by comparing the exper- 
imental ARPES result 11 (upper panel) with the theoretical 
one (lower panel), which has t — — 0.3i and molecular rota- 
tion angle 9 — 64°. The energy unit has been converted to 
eV using t = 2.72eV. The momentum unit is a 1 . 
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FIG. 6: (Color online) Comparison of DOS from QMCPT cal- 
culations and experimental PES data^i for C60 (upper panel) 
and K3C60 (lower panel) hexagonal monolayer. The QMCPT 
energy scale has been converted to electron volts by setting 
t = 2.72eV. 
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FIG. 7: Band dispersions along TM direction with different rotation angles 6 for a 2D hexagonal K3C60 superlattice from 
QMCPT calculations. The other parameters are fixed: U — At, t — —O.St, and T = O.li. The momentum unit is A -1 . 
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Fermi energy for the K3C60 monolayer. In addition, we 
find small peaks at around — 2eV and valleys at around 
— 4eV for both theoretical and experimental PES curves. 
Thus, although we do not claim that U = 4t corresponds 
to a fit to the data, this value of U is consistent with the 
PES data. Smaller or larger values of U would lead to 
smaller or larger values of the gap for neutral Ceo mono- 
layers. Unfortunately, to go beyond this by varying the 
value of U would, at present, require a prohibitively large 
amount of computer time. 

We comment that the above-determined value of t 
is consistent with earlier TB and LDA calculations 23 , 
where Satpathy et al. found good agreement between 
TB and LDA band structure calculations for a fee Ceo 
lattice. The magnitudes of their TB hopping integrals 
between two nearest C atoms on two NN Ceo molecules 
are, depending on the relative orientations of NN Ceo 
molecules, in the range of (0.50eV, 0.82eV). From our 
calculation, the maximum C-C hopping integral between 
two NN Ceo molecules is about 0.51eV (calculated from 
t Oin,Jj)- Therefore, our estimated TB hopping integrals 
agree with Satpathy et a/.'s calculations. 

In our formalism, it is easy to explore the variation of 
the monolayer band structure with different Cgo rotation 
angles 9. In fact it was necessary to do just this to de- 
termine the orientation which best fit the ARPES band 
structure. Fixing t = —0.3t, we carry out QMCPT for 
Ceo molecular rotation angles 9 = 0, 10°, • • • , 70° for a 
K3C60 monolayer. The band dispersions along the TM 
direction are shown in Fig. [7] We see that the band struc- 
ture is quite sensitive to the molecular rotation angle 9. 
For example, at 9 = 30°, the band is inverted with a 
band minimum at T and a maximum at M. At 9 = 40°, 
the band becomes almost insulating. These are quite in- 
teresting, and might be verified by experiments if a way 
can be found to control the molecular orientation. 



IV. CONCLUSION 

To summarize we have calculated theoretical single- 
particle excitation spectra for single Cqo and CgQ 
molecules as well as for hexagonal monolayers of Cgo and 



K^Ceo- The calculations were carried out by the QM- 
CPT technique ; 16 ' 17 ! 18 using a Hubbard model for each 
molecule, and simple electron hopping between molecules 
of the 2D Ceo superlattice. This new computational 
method is particularly well-suited to composite systems 
of fullerenes. Using a plausible value, U = 4t, for the 
on-site Coulomb interaction U, we fit the effective NN 
molecular hopping amplitude, t , and the Ceo molecu- 
lar rotation angle 9 to ARPES and XPD data. We find 
that the electronic spectra are extremely sensitive to the 
molecular orientation angle 9. With these fitted param- 
eters, we find good agreement with photoemission data 
for the electronic density of states. The C-C hopping in- 
tegrals between NN molecules obtained from the fits are 
consistent with earlier TB and LDA calculations. 23 

The significance of the present calculations is to show 
that the electronic spectral data obtained so far are com- 
patible with the relatively simple model of electronic cor- 
relations contained in the intramolecular Hubbard model 
with simple interball hopping, in much the same way as 
they have also been shown by other o 11 ' 12 to be compat- 
ible with LDA/phonon calculations. Further measure- 
ments and calculations will be required in order to refine 
our understanding of the dominant correlations in these 
systems. 
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